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For five- dimensional braneworlds with an Si/7j2 orbifold topology for the extra dimension x^ , we 
discuss the validity of recent claims that a gauge exists where the two boundary branes lie at fixed 
positions and the metric satisfies g^s = dsgss — where fi labels the transverse dimensions. We 
focus on models where the bulk is empty apart from a negative cosmological constant, which, in 
the case of cosmological symmetry, implies the existence of a static frame with Schwarzschild-AdS 
geometry. Considering the background case with the branes moving apart after a collision, we show 
that such a gauge can be constructed perturbatively, expanding in either the time after collision or 
the brane velocity. Finally we examine how cosmological perturbations can be accommodated in 
such a gauge. 



I. INTRODUCTION 

The large number of recent publications on braneworld 
cosmology is a testament to the richness and diffi- 
culty of this comparatively new field (see, for example, 
[H Q> H> Hi Hi ly|)' Despite the large gap in complex- 
ity between possible String/M-Theory realisations and 
the greatly simplified models currently being investi- 
gated, even the simplest models are extremely rich in 
new ph ysic s. Ma ny s uch five-dimensional toy models 
0, S MM m El [il H are motivated by the dimen- 
sional reduction of Hofava-Witten M-theory >T5l| to five 
dimensions but miss out most of the resulting large num- 
ber of bulk degrees of freedom. The extra dimension is 
assumed to have the topology of an 81/1^2 orbifold, with 
two 3-branes residing at the boundaries of the spacetime 
at the two fixed points. 

The increase in complexity from a purely four- 
dimensional Universe to one confined to a brane embed- 
ded in a five-dimensional bulk is considerable. We shall 
therefore mainly focus on the Randall-Sundrum model 
Q, where the bulk is empty apart from a cosmological 
constant. We assume that the branes contain no matter 
and carry equal and opposite tensions finely tuned with 
respect to the bulk cosmological constant such that their 
induced cosmological constants vanish, although the con- 
clusions of this paper will not depend sensitively on these 
assumptions. In the background case, where the trans- 
verse dimensions are assumed to possess cosmological 
symmetry, the lack of non-gravitational degrees of free- 
dom in the bulk then implies the existence of a Birkhoff 
frame in which the bulk is manifestly static with Anti-de- 
Sitter(AdS)-Schwarzschild geometry, parameterised by C 
(related to the black hole mass) and the lengthscale L. 
The velocities of the branes are then determined entirely 
by their positions, and the system can be solved exactly 
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even in the presence of matter on the branes; the only 
deviation from four-dimensional Einstein gravity is the 
presence in the Friedmann equation of quadratic terms 
in the brane stress-energy and a dark radiation term pro- 
portional to C, which contains all of the gravitational 
effects of the bulk. 

When one comes to include cosmological perturba- 
tions |l(i I13, however, the problem rapidly becomes 
intractable, largely due to the considerable increase in 
number of the gravitational degrees of freedom in going 
from four to five dimensions. To achieve any sort of an- 
alytic result on the propagation of braneworld perturba- 
tions one needs to resort to approximations. The simplest 
is to work to lowest order in the brane velocities, assum- 
ing also that the energy content of the branes is much less 
than their tensions. The result is a four-dimensional ef- 
fective scalar-tensor theory li,'!!,'!!!!!!!!!!!!!!!!!, 
valid at low energies, formulated for example in terms of 
the induced metric of one of the branes. The scalar de- 
gree of freedom, known as the radion, corresponds to the 
distance between the two branes. 

Attempting to go beyond the very limiting low-energy 
approximation requires a method of solving the five- 
dimensional Einstein equations. While a great deal of 
progress has been made recently [H IH Hi IM El S 
iSl III HI HI Hi, the problem is still very difficult. The 
five-dimensional analysis is rendered simpler by choosing 
a gauge in which the branes are located at fixed positions. 
The further assumptions that the coordinate x^ param- 
eterising the extra dimension is chosen such that both 
the cross terms 55^ (where ^ = 0, f , 2, 3 labels the trans- 
verse time and space dimensions) vanish and 95(755 = 
result in a particularly useful gauge, henceforth referred 
to as a simple fixed-brane (SFB) gauge. Several publi- 
cations m I33, in, I33 have already assumed, without 
proof, the existence of SFB gauges to derive effective 
theories with different, more useful regimes of validity 
than the low-energy scalar-tensor theory. In particular 
they focused on the small-radion limit pertinent to re- 
cent attempts to model the Big Bang as a brane collision 
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This paper is organised as follows. In S|n]we discuss 
the general problem of constructing SFB gauges, special- 
ising in ijllll to the special case of background spacetimes 
with cosmological symmetry. For the case of two branes 
moving apart after a collision, initially considered in the 
Birkhoff frame, we construct in ^IVI the required coordi- 
nate transformation as a power series in the SFB time t, 
equivalent to an expansion in the radion. We show that 
such a gauge can be found correct to arbitrary orders 
in t, although there is in general no analytic expression 
for the terms in the series and one must resort to nu- 
merical methods. We also comment on the convergence 
of this series and its relation to the horizon structure of 
the bulk. A second method is derived in ^V\ namely an 
expansion in \/C, which is equivalent to the velocity of 
the branes. This turns out to be simpler, with analytic 
expressions available at all orders. Finally we discuss cos- 
mological perturbations in WVI\ and show the existence of 
the gauge transformations which will accommodate gen- 
eral perturbations into an SFB gauge. 



II. GENERAL FORMALISM 

We take our spacetime to be an 81/^2 orbifold with 
two boundary branes at the fixed points. We assume 
first of all the existence of a Gaussian- Normal coordinate 
system relative to the positive tension brane at y = 0, 
with metric 



ds^ = dy^ + h^^ {x, y) di^'di''. 



(1) 



Throughout this paper Greek indices will run over 
0,1,2,3 and denote the transverse spatial x' and tem- 
poral x" coordinates, and a, 6, c etc. are five-dimensional 
indices running over 0,1,2,3 and 5 for the orbifold di- 
mension. The location of the positive-tension brane at 
y+ = is not crucial for the discussion, of more impor- 
tance is the form of the metric with 355 = 1 and (75^ = 
for the sake of simplicity in the following argument. 

Though the coordinate system |^ is particularly sim- 
ple, the second brane will, in general, lie on a spacetime- 
dependent trajectory y — y_(x). Implementing the Israel 
junction conditions on a moving brane is a complicated 
procedure 44J : far simpler is to use a coordinate system 
in which the branes lie at fixed positions, at the cost of 
losing the simple form of the metric ^. A general coordi- 
nate transformation to bring the branes to fixed positions 
will introduce off-diagonal terms 5^5 and a non-trivial 555 
to the metric. However, the aim of this paper is to show 
the existence of a 'simplest possible' SFB gauge where 
both the non-diagonal metric terms vanish and 355 has 
no dependence on the extra dimension y 



ds^ = d{x)^dy^ + g^,^{x,y)dx''dx'' , 



(2) 



a four-dimensional scalar field which plays the role of the 
radion, giving the proper distance between the branes on 
curves of constant x'^. Note that, as discussed in [33|, this 
definition of the radion is frame dependent. As we shall 
see in ijlll Al this function is expected to be monotoni- 
cally increasing (for the homogeneous cosmological back- 
ground) for the case of two branes moving apart after a 
collision, and vanishes only at the surface of collision. 
Defining the diffcomorphism 

x'^^r{x'^,y),y^g{x'',y), (3) 

we must then find C°° functions /^ and g such that 

do df 

-^d^g + K,{r,g)-L-d^r - 



(4) 



dg 



a/M df^ 



in which the branes lie at the fixed positions y = and 
y = 1. The metric coefficient d{x) can then be viewed as 



gj +^..(r,5)^^ = 4^"f>o. (5) 



The new metric is then of the form (0). .9(2^) niust 
then satisfy boundary conditions related to the positions 
y±{x) of the branes in the original frame (^ 

g {x, 0) = y+ {fix, 0)) , g {x, 1) = y_ {f{x, 1)) (6) 

in order to bring the positions of the branes in the new 
frame to y = for the positive-tension brane and y = 1 
for the negative-tension brane. Note that there is a resid- 
ual gauge freedom of a;-reparameterisation which leaves 
both the form of the metric and the brane positions in- 
variant. 

This system of non-linear, partial differential equations 
is extremely difficult to tackle in general. However, since 
we are ultimately interested in cosmological braneworlds, 
we will at first only need to consider these equations for 
the homogeneous case, effectively reducing the dimen- 
sionality of the problem from five to two. Furthermore, 
this paper is only concerned with Randall-Sundrum style 
braneworlds where the bulk contains only a (negative) 
cosmological constant; staticity of the bulk then follows 
in the homogeneous case from Birkhoff's theorem. 

This static Birkhoff Frame is of the form (^ and we 
investigate its properties in the next section. We will de- 
scribe two methods of constructing the coordinate trans- 
formation satisfying Q and (|SJ| for the background, and 
then turn our attention to cosmological perturbations. 



III. BACKGROUND FORMALISM 

A. Birkhoff Frame Dynamics 

In this section we re-derive some standard results on 
the propagation of branes in the static Birkhoff frame 
background 33, 34, 45], with metric 



ds' 



dY^ - n^{Y)dT'^ + a^{Y)dy? 



(7) 



where 



n\Y) 



,-2Y/L _ 



C 



3^^/^ (8) 

4 
a'{Y) = -n{Y)/L (9) 



The bulk geometry is Schwarzschild-AdS, with the AdS 
lengthscale L related to the cosmological constant A by 
A = — -p- and the parameter C proportional to the mass 
of the bulk black hole. We assume that the brane tensions 
are fine-tuned to the usual values A± = ±-|; where the 
five-dimensional gravitational constant has been set to 
unity for convenience. The spatial curvature is assumed 
to vanish for simplicity. 

In this frame the branes will not be static, but fol- 
low trajectories Y±{T). Assuming for the moment that 
the branes are empty, the Israel junction conditions are 
simply 



K\ 



a± 






(10) 



where K^ are the extrinsic curvatures of the branes. 
The space-space components of this equation then give 
the first-order equations of motion [SJ, yj] 



dY±Y 
dT 



C- 



(11) 



where a± — a{Y±) etc. Since we are interested in brane 
collisions we assume that dY±/dT are of opposite signs, 
specifically that dy+/dT < and dr_/dT > 0. The 
proper time t± on the branes can then be determined, 
giving the world velocities and accelerations as 



dY± 
dt± 

D^Y± 



C 



a±n± 
a± 



n±L 



d7± 
dt± 



D4 



a± 



C 



n^a±L 



(12) 
(13) 



The above equations imply that, for empty, fine-tuned, 
spatially-flat branes, one cannot achieve a consistent em- 
bedding into the negative mass Schwarzschild-AdS space- 
time with C < 0. It also implies that, under the same 
conditions, the C = solution is the trivial case where 
the branes remain static and the frame reduces to the 
usual AdS warp-factor profile, although from (|13|l we see 
that the branes are still accelerating. 

We shall therefore henceforth only consider the case of 
positive spacetime mass C > 0. In this case there is a 
horizon in the bulk where n(Y) vanishes, given by 



from this coordinate system only to be valid up to the 
point where one of the branes hits the horizon (in other 
words, for the physical slice of spacetime bounded by the 
two branes to develop a horizon). Since T-translation 
is a symmetry of the metric (|7|l we can assume that 
the brane collision occurs at Y = Yq, T = 0. Fig. 1 
shows the branes' trajectories in the Birkhoff frame with 
L = C = 1, lo = 0. Note that a brane crossing the hori- 
zon will do so in finite proper time but will require an 
infinite AT, since from Hll|l we see that the brane ve- 
locity with respect to the Birkhoff time T tends to zero 
as the horizon is approached. The corresponding picture 
as a Penrose diagram is illustrated in Fig. 2, showing the 
region of the extended Schwarzschild-AdS spacetime cor- 
responding to the coordinate patch (Q). In this picture, 
the horizon at Y = l^or can be interpreted as the event 
horizon of the bulk black hole. Note that the physical 
spacetime is just the slice bounded by the two branes 
(or, rather, two copies of it due to the orbifold nature of 
the complete spacetime). 

Importantly, it follows from (|ll|l and l|12|) that the 
Friedmann equation on the brane is simply Hj_ = 
C/jL'^a X). i .e. the background geometry is finite at the 
collision [43|. The event is, however, singular from a five- 
dimensional perspective, in that the fifth dimension mo- 
mentarily vanishes. In general perturbations will diverge 
due to the blue-shifting of their energies as the radius of 
the fifth dimension vanishes, although several authors re- 
cently have claimed that this need not be an obstacle to 
the propagation either of cosmological perturbations |4ll | 
or, more fundamentally, string excitations |40ll46| across 
the bounce. 



B. Gauge Transformation 



We now simplify the analysis of 
neous, static background of illll Al 
ponents iQ and the brane trajectories 



l] to the homoge- 

The metric com- 

((TT)l do not. 



by homogeneity, depend on the spatial coordinates x^. 
We can therefore just leave these coordinates alone, 
/* {x, y) = x^ , and focus on the two-dimensional Y, T por- 
tion of the metric. As before, we consider a diffeomor- 
phism 



T = /(y,i), Y^giy,t). 



(15) 



The conditions l@J and ((SJ on the form of the new metric 
become 



9v9t = nig) fyft 
9y 



' = n{gffl + d{t)\ 



(16) 
(17) 






(14) 



The coordinate system ((JJ will only be valid in a par- 
ticular coordinate patch bounded by the horizon. One 
would therefore expect any results in this paper derived 



givmg 

ds^ = d{tfdy^ - b {y, tf dt^ + a {y, tf dx^ 
where 

h{y,tf^n{gfj^^gl 



(18) 




0.4 



FIG. 1: Background brane trajectories for L = C — l,Yo — 0. 
The trajectory of the positive- (negative-) tension brane is 
the solid (dashed) hne, and the event horizon at Y = Yhor is 
denoted by the solid vertical line. 




FIG. 2: The Penrose diagram corresponding to Fig. 1, i.e. for 
the case C > 0. The shaded area denotes the region of ex- 
tended Schwarzschild-AdS covered by the coordinate system 
Q. The trajectory of the positive- (negative- ) tension brane 
is denoted by the solid (dashed) line. The physical region 
of spacetime is that bounded by the two branes. The diago- 
nal lines bounding the coordinate patch represent the horizon 
Y — Yhoi- Region I therefore represents the portion of the 
future of the collision which the Y, T coordinates cover. 



It is at this point that the assumption that the metric 
is static (i.e. that we are concerned only with Randall- 
Sundrum-style braneworlds, with no non-gravitational 
degrees of freedom in the bulk) comes in. This implies 
that n is only a function of g, rather than n(/, g) as would 
be expected more generally. One can then eliminate / 
entirely using (|15|l and H17|) . giving (up to an arbitrary 



sign) 



ft 



1 



nig) 
1 



9l 



d{tf 



9y9t 



<9)^g2-dity 



(19) 
(20) 



The Jacobian of this coordinate transformation must not 
vanish, which implies from H19|) and (|20|l that 

d{t)^gt y^ 0. 

We can therefore conclude that gt ^ away from the 
collision which, by convention, we shall assume occurs at 
T = 0. In the (y, t) plane, the collision corresponds to 
the vanishing of the metric component d{t), and also to 
the curve f{y,t) — 0. This curve has zero proper length, 
lying at constant t, and wlog we can take the collision to 
correspond to the line i = on which / vanishes and g 
takes the constant value Iq- 

DIT31 



As we saw in illll Al the distance between the two 
branes is a monotonically-increasing function after the 
collision. This implies that one can exploit the t- 
reparameterisation of the metric to fix 



d{t) = t. 



(21) 



It turns out that this gauge choice ensures the existence 
of a perturbative expansion of the functions / and g in 
powers of t, as well as corresponding to the linear ap- 
proximation for close branes used in [3^. 
Continuity implies that fyt = fty, which from (|19|) and 
(Pn|l reduces to a second-order partial differential equa- 
tion for g [y, t): 



9l , n'{g) 



9t9vy + t ^ , ^^^^ y^.yy 



9t9l - t^9t) = 



(22) 



Writing (HJl) as 



dT 



Tv{Y±{T)), 



we find 



d 



Y+{fiOA) = -v{g)ft\, 



y=a 



9t\y=0 - Q^^ 

9t\y=, = §-Y-if{l,t) = + vig)My^, 
Substituting l(^ for ft we find, for example. 



(23) 
(24) 



9t\ 



v{g) gygt 



y=o 



<9)^l~rT^ 

n{gft^ 



y=o 



y=o 



■^y\y=o n{gf-v{gf 

Using lO and (|11|) then gives the boundary conditions 

aig) 



5yly=o,i 



= i 



n{g) 



(25) 



y=o,i 



together with 



B. Second order 



g{y,t = 0)=Yo. 

Note that the sign of gy in (|25|) is specified by the choice 
ft > 0, gy > 0, i.e. that the orientation of the new coor- 
dinate system is the same as the old one. The fact that 
the branes are moving in opposite directions, not imme- 
diately apparent from the boundary conditions since the 
— sign in 123|l has been lost, follows from the choice of 
gauge lEIJ. 



IV. PERTURBATIVE EXPANSION IN t 

We now proceed to solve (I22|l subject to the boundary 
conditions H25|l as a perturbative expansion in powers of 
t. Specifically, we write 



g(y,i)=ro + ^ff„(2/)r 



(26) 



where the boundary conditions H25() will specify the 
derivatives of g^ at y = 0, 1 in terms of the values there of 
the lower-order gt [i < n). (|22|l becomes a set of ordinary 
differential equations for the 5„; the aim therefore is to 
identify the correct initial value 5n(0) for which the two 
boundary conditions are simultaneously satisfied. Once 
g{y,t) has been determined to the required order, f{y,t) 
can in principle be recovered from (|20|l with boundary 
condition /(y, 0) = 0. 



A. First order 



The 0{t) contribution to (j^ gives 



9i9i - .gf + 1 = 



with l|25(l implying 

5;(o)=5;(i) = - 

where ao = a(yo) etc. We then obtain the result 

5i(y) = T^sinh(a(2?;- 1)), cosha= — 
Za no 

The corresponding result for / is 

1 



fiiy) 



2anQ 



cosh (a {2y — 1)) , 



(27) 



(28) 



(29) 



(30) 



This is just the relevant coordinate transformation in the 
approximation of constant brane velocities, as detailed in 



Before examining the existence of gn for arbitrary n, 
we consider the simpler case of n = 2. 



5152 - '^9192 + 2525" 



n'jYo) 
n{Yo) 



{919? 



91) =0 (31) 



Expanding out 

K{g{y,t))^K{Yo)+tg,{y)K'{Yo) + 0{t^) 



(32) 



for an arbitrary analytic function K, we identify the 
boundary conditions as 



givmg 



5^(5 = 0,1)^ (^)' 



52(1) - -52(0) 



.9i(y = o,i) 



Y=Yo 



C sinh a 
Lan^ttQ 



(33) 



There is no simple analytic general solution to this equa- 
tion, so one must proceed numerically. The task is to 
find the value of ^2(0) such that, after integrating H31|) 
subject to the initial conditions at y = 0, 52(1) is given 
by H33() . for example using a bisection procedure. The 
existence of g2 is discussed below. 



C. General n 

For tuned tensions, i.e. the boundary conditions l|25|l . 
the gn all have definite parity (—1)" about y = 1/2. To 
show this, we substitute the power series ()26|l into l|22|l 
and consider the coefficient of i", giving 



515" 



'^9Wn 



ng'i9n 



Pniv) 



(34) 



where P„ (y) is composed of a finite sum of products of 
the gi and their first or second derivatives, but only con- 
taining gi for i < n (c.f. H31() ). Examining the contri- 
butions from the various terms in H22|l individually, we 
write symbolically (i.e. ignoring numerical constants) 



9t9yy 



fe=0 




5fc+i- 




where the symbol D denotes the terms in the coefficient 
of i" which do not contain gn, i.e. the contribution to 
the RHS of (|5Hl . Working under the inductive hypothesis 
that gm has parity (—1)™ for m < n (which is obviously 
true for n = and n = 1), we see that each summand 
has the parity (—1)"+-'^. Similarly, 
tribution 



-f^gy makes a con- 



1 



n-l 

^dY gin 9'. 



•Jy 



m=2 



n+1 — m 



Again, the parity of these terms is (—1)"'+-'^ by the induc- 
tive hypothesis. Next, we examine gtg^ — t^gt, which is 
multipHed by a function of g in 122|l , so we will need all 
orders up to n in its expansion. Wc find that 

oo / k—1 k — m \ 

5t5y ^ X!^'' ( X! X! 9m g'p g'k+1-rn-p j 
fc=2 \m=l p=l / 

n 

fe=2 

where Hk{y) is a function of parity (—1)''+^ involving 
only gm for m < k; there is no contribution at order i" 
or below from g^- t^gt is obviously also of this form, so 
we just need to look at the n'{g)/n{g) term. Any func- 
tion K{g) with a power series expansion can be written 
symbolically in the form 



fc=l ym=l ij>O.Y,ij=k 

(35) 
which can be re-expressed as 

K{g) ^ (J2t^Gu{y)\ +gn{y)e + 0{e+^) 

where Gk{y) is a sum of products of gi with i < k, with 
total parity (—1)''. Putting this all together, we see that 
the RHS of 134|l can be written as 

n— 1 n— 1 

"n '^ 2_^ 9m9n+l-m + 2_^ 9m 9n+l-ni 



y^ Gjn{y)Hn-„i{y)- 



(36) 



m=0 



Each summand in each term has the same parity, 
(—1)"+^. Therefore, looking back to IJ23J, we see that 
the whole equation is invariant under parity if gn has 
parity (— 1)". Whether or not this is the case depends 
on the boundary conditions; in the fine-tuned case, the 
boundary conditions are just that 



?n(y*) ^ Fn~i{y*) 



where y* 



0, 1 and 

a{9{y,t)) 
n{g{y,t}) 



E Pr.{yr 



Note that this equation implies that 5^ at the boundaries 
is specified entirely in terms of the values there of the gm 
for m < n. From l|35|) . Fn-i{y) has parity (— 1)"~^. 
Hence 



and the inductive hypothesis is proven. One can now 
show the existence of <?„ for n odd - since the function 
will be odd, it and its second derivative must vanish at 
y = 1/2. The equation of motion will then prescribe the 
value of g'^{l/2) which, together with 5„(l/2) — 0, will 
give a unique solution. 

For even n, one must use a more general argument. 
We are searching for the value of 5n(0) which gives a 
particular value for 3^(1). Taking (7„(0) = A^ for A^ very 
large, but 5^(0) to its finite, prescribed value, the LHS 
of (|34|) will dominate over the RHS apart from at the 
finite number of points where either gn or its derivatives 
vanish. Therefore, 



gn{y)^NK{y) as N ^ g^i^) 



-CX) 



(38) 



g'nW = i-ir-'g'ni^) 



(37) 



where /i„ is the solution to the initial value problem 
9ifn - 29'Jn + ng'Ifn = 0, /„(0) - 1, /;(0) = 

Provided that /i^(l) 7^ (which can be shown numeri- 
cally for specific values of n), we then have that (7^(1) -^ 
±00 as (?„(0) -^ +00, g'n{l) -^ Too as 5,1 (0) -^ —00. 
Hence, by continuity, there will be (at least one) inter- 
mediate value of gn(0) where g^(l) takes the required 
value. Hence we can always find a gn satisfying both 
(|34|l and the two boundary conditions for 5^ at y = 0,1 
following from H25|) . 

The above argument is not sensitive to the precise na- 
ture of the boundary conditions; in other words, in the 
case of general brane tensions and brane matter content, 
the boundary conditions (|25|l would be modified, but one 
would still expect to be able to find suitable functions 
gn{y) to perform the coordinate transformation to the 
SFB frame. 

However, the entire analysis is dependent on the ex- 
istence of a static bulk frame such as Q (though the 
precise nature of the functions a and n will, again, not 
affect the argument). It is not clear if this result could 
be generalised to more complicated braneworlds where 
the bulk is not empty, such as dilatonic or su pergravity - 
inspired models with bulk scalar fields pTllTlfll llll47| . 



D. Numerical results and convergence 

The horizon structure of the bulk has ramifications for 
the regions of validity of the coordinate system and the 
radius of convergence of the power series (PH|l . For C > 
there is a horizon, whilst for C < there is a naked singu- 
larity and no consistent embedding of flat, empty, finely- 
tuned branes. The C = case, for fine-tuned branes, 
represents the trivial evolution where the branes remain 
fixed at their original positions, and the Birkhoff frame 
automatically satisfies the required gauge conditions |(2J). 
We therefore only need consider the case C > 0. 

The first six non-trivial (/„ have been computed nu- 
merically using a bisection procedure to search for (?n(0) 
in the specific case oiYo~Q, L = C = 1. The resulting 
52 to (7e are plotted in Fig. 3. 
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FIG. 3: g2 to ga plotted for L 
alternating parities. 



1,C = l,yo = 0. Note the 
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FIG. 4: Comparison of the brane trajectories computed us- 
ing the coordinate transformation function g{y, t) up to 0{t^) 
with the exact result. The dotted curves correspond to the 
exact trajectories as in Fig. 1, showing highly accurate agree- 
ment with the numerical results. For comparison the corre- 
sponding trajectories computed working only to first order in 
t are also shown for the same range of t, with clear differences. 



Fig. 4 gives the calculated brane trajectories back in 
the iX:T) frame; for example, 5(0, t) gives Y^{t), and 
integrating 1(201 with d{t) = t gives /(O, t) = T_ (t), hence 
one can plot the trajectories Y±{t) and compare them 
with the exact solutions. The accuracy of the agreement 
is high, as can be seen from the Fig. 5. Also plotted 
are the trajectories calculated using only the first-order 
expression for g(y,t), with the same range of t; these 
clearly differ wildly from the higher-order results as the 
horizon is approached. 

Depending on the initial conditions, either the positive- 
or negative-tension brane will hit this horizon at a finite 
value of t = th- For the initial conditions chosen above, 



this corresponds to th « 0.36. At this point, n{g) will 
vanish at the relevant boundary, and so (|22l) will break 
down and one would not expect the coordinate system 
(|18|l to be valid beyond this value of t. In particular, one 
would expect this value of t to correspond to the radius of 
convergence of the series (|26|l . It is clear from Fig. 3 that 
the boundary values of gn{y) are increasing in modulus 
as n increases, which is a manifestation of the lack of 
convergence of the series for large t in this case. 

In the fine-tuned, empty-brane case with boundary 
conditions H25|l . the symmetry of the problem (i.e. the 
manifest reflection symmetry under y ^ — y in the tra- 
jectory plots) implies that the same coordinate system 
will also be valid before the collision, for —th > t > 0. 

The advantage of an expansion in the SFB time t is 
that it is ideally adapted to working in the limit of close 
branes, a regime of considerable interest recently in the 
literature pE S ii ^ 48, 49 J . In effect, an expansion 
in t is equivalent to an expansion in the distance between 
the two branes. Working even just to first order in i or d 
whilst making no other approximations, one can obtain 
interesting, non-perturbative results unattainable using 
the moduli space approximation |33 . yjf. Whether or 
not the full series H26|) converges is then not important if 
one simply wants to proceed perturbatively in t. 



PERTURBATIVE EXPANSION IN C 



An alternative method of finding the required gauge 
^ is to consider an expansion in the velocity of the 
branes l|ll|l . Specifically, we define the parameter 



h = ^/C 



(39) 



and expand in this rather than in t. This amounts to 
an expansion about the C = AdS spacetime where the 
(finely-tuned) branes are at fixed positions in the Birkhoff 
frame Q. This method turns out to be conceptually 
and analytically simpler than the expansion in t, though 
it is of less direct use for the near-brane regime, being 
an expansion about a static configuration. In this sense, 
working to linear order in h will reproduce the results of 
the moduli space approximation. 

As discussed below, we will in fact expand about the 
state where the branes are static and coincident, in order 
to simplify the analysis somewhat and to tie in with the 
rest of the paper assuming that the branes are moving 
apart from an initial collision. It would not, however, 
be difficult to generalise the results to a different initial 
configuration. 



A. Formalism 

Firstly we define a new coordinate Z — e~^'^ , in terms 
of which the BirkhofF frame metric can be written 

ds^ = -^dZ^-n{ZfdT^ + a{Zfdx^ (40) 



n{Zf = a{Zf- 



a{Zf 



.(+) 



and 



For ft, = the branes are static, a,t Z 
Z = Zq'. Since we are interested in modeUing the af- 
termath of a coUision we shall take Zq = Zq. Whilst 
simplifying the resulting analysis considerably one must 
then remember that, for /i = 0, the metric will be sin- 
gular. Therefore we shall expect gyy in the SFB frame 
(|18|l to vanish for h = 0. Bearing this in mind, we then 
perform a coordinate transformation as before, keeping 
the i-gauge general for the time being (as opposed to the 
choice d{t) = t of the previous section). Expanding in h, 
we write 

Z{y,t) = Zoh+f;ft"Z„(j/,i)j (42) 

oo 

T{y,t) = t+^ft"r„(y,t), (43) 

n=l 

where for h — we can set T — t since the branes are 
already at fixed positions. We then solve 

dy9yy = 9yt = 

order by order to determine the functions Zn{y,t) and 
Tn{y,t), together with the boundary conditions (|11() . 



d^± ^ „ h n{Z±) 
dT ^ ^La{Z±)^ 



B. First order 



(44) 



In the (y, t) frame, the metric becomes 



2 _ I ^ ry2 



d^' = ["^Zy-niZrT^jdy' 

-2 { ^ZyZt - niZfTyTt ) dydi 



(45) 



ill 

Z2" 



2ti2 1 1.2 I „/'7^2J„2 



Zf - n{ZyTf dr + a{Zydx- 



and so, using (|42|l and (|43|) . we obtain 

dygyy = 2 (^L^ Z I Z '^ - Z^ T^T'^^ h^ + (h^) (46) 



gyt = -Z^tih + Oih^ 



Next we examine the boundary conditions (|44|l . 
dZ+ _ dtZ{0,t) 



dT 



dtT(0,t) 



= ZoZi{0,t)h + 0{h^) 



4^4-«f^' 



givmg 



(■41) and, similarly. 



Zi(0,i) = 



Zi(l,t) 



1 



LZo 
1 

Tz^- 



(48) 



(49) 



(|47|l gives the simple condition that Ti = and, choosing 
the collision to occur at t = 0, implies 

T,iy,t) = 0. (50) 

The condition ipi)| on gyy reduces to 

Zj' = 

which, together with the boundary conditions H48|l . H49|) 
and ^i(y,0) = Zq, gives 

Ziiy,t) 



LZo 



(1 - 2y) 



(51) 



Therefore, to 0(h), the required coordinate transforma- 
tion is 

Z{y,t) = Zo + ^{l-2y) + 0{h') (52) 

Tiy,t) = t + Oih"). 

The metric (|45f) is then given by 
/ f2/,2 \ 



ds = 4— 2- + 0(ft-') d2/"+0(ft")dydt 



72 
^0 



.Znt 



+ {ZI + 2-^ (1 - 2y) ft + 0(ft2) j (-dt^ + dx^) 

which, to 0(ft), is of the required form |(5J. As a check, 
we can calculate the proper times on the two branes, 



di± = ( Zo ± -^ + 0(ft2) ) dt 



giving the Hubble parameters as 

ft 



H. 



± 



LZl 



Oih^). 



(54) 



(55) 



The function d{t) , giving the proper distance between the 
branes along lines of constant y, can be read off from H53(l 
as 



which gives 



9f 

Za 



H.-±^^^m^oih% 



(56) 



(57) 



(47) in agreement with the low-energy result of |3j 



C. Higher orders 

The analysis at higher orders is straightforward. The 
constraint dyQyy = reduces to an equation of the form 



fin 



^"(y,0 = E4r^^ 



(58) 



where the superscript means that the sum is finite, i.e. a 
binomial in y and t. This can easily be integrated to ob- 
tain Zn{y, t) up to two arbitrary functions of t, which are 
fixed by the boundary conditions H44|) and the require- 
ment that Z{y,0) = Zq. Similarly, the other constraint 
gyt = is of the form 



K{y,t)^Y.€'t'v 



(59) 



J J 



giving Tn{y, t) up to an arbitrary function of t. This is a 
manifestation of the i-reparameterisation invariance, and 
can be fixed for example by setting T{t, 0) = t. In this 
particular gauge, the next two sets of functions are 



Z2{y,t) 



2i2 



L'Zi 



y{y- 1) 



2t 



T2iy,t) = --^y(l-y 



Z3iy,t) 



T3iy,t) 



(60) 
(61) 







-2 I .2 r^2\ I y y 



,3^.^4(X^+.%^)^^-^j (62) 



2t'y 

'dLZl 



(9-12y + 8y2 



(63) 



This transformation satisfies all the required boundary 
conditions, as well as the conditions on the metric ^, 
up to 0{h^): 



9yy - 4^ + ^^^ 

4\ 



It^Z^- 19L^) +0{h^) 



9yt = 0{h^) 



In order to prove the relations l|^ and (|CT|) . we use in- 
duction. Specifically, we aim to prove that the general 
form of Zn and T„ is a binomial function of t and y of 
the form 



Z^{y,t)^tY,^t!y't' forn>l 



(64) 



i,j=0 
fin 



T„(2/, i) = y t E b[fyH^ for n>2. (65) 



From (|51|l and (|60I63|I . we see that this is indeed the 
case for n < 3. We now assume that we have found the 



Zfn and Tm for m < n, satisfying the required boundary 
conditions. We then look at the constraint that gyy is 
independent of y, which becomes 



ZZyZyy = Z ^ + — {n{Z)n' (Z) Z^ Z yT^ + niZ)Z%Tyy) . 

(66) 
We aim for an equation of the form H58|) and, examining 
the LHS of this equation, we see that the lowest order in 
h at which Z" appears is at 0(/i"+^), 

ZZyZyy D {z^ Z{ Z'^ + t^ B {y , t)) h''+^ 

^Z^oZ:: + t'Biy,t)y-+' (67) 

where the symbolic notation B{y,t) represents a bino- 
mial function in t and y depending only on the binomial 
functions Z^ for m < n. The factor of t^ multiplying 
this function is displayed explicitly since we need to keep 
track of powers of t to ensure the RHS of H58|) is indeed a 
binomial with no negative powers of t. Now starting on 
the RHS of ^, we note that 

OO 

Z3 = Z3 J2 h'^^'+^Z'^ZlZ'^ 

a,6,c— 1 

and hence, at 0(ft,"+^), does not contain Z^, only the 
lower Zm- Hence we can write symbolically 

Zl^t^B{y,t). 

Both Ty and TyTyy have a factor t^ h* in front, and can 
be written with the inductive hypothesis as 

T^, TyTyy ^ t^ h^ B {y , t) + t^ O {h''+^) . (68) 

Furthermore, as in l|35|) . a general analytic function of Z 
can be written 



f(z) = f(z,)+Y.h^\Y^ E 



(69) 



fe=i 



, m=l 'ij>0,J3 ii=k 



Combining (|68l) and (|69() , we see that the term with coef- 
ficient l/L^ on the RHS of ^ contributes, at 0(/i"+i), 
a term of the form PB{y, t). 

Hence, putting all of the above together, we have 
shown that Z„ is given in terms of the Zm for m < n 
by an equation of the form (|58ll with an additional pre- 
factor of t. This equation can always be integrated to 
give 



Zn {y,t) = tBn {y,t) + an{t)y + bn (t) 



(70) 



where Bn{y, t) is a binomial expression depending on the 
Z,„ for m < n, and a„(i) and bn{t) are arbitrary func- 
tions, which are then uniquely fixed by the boundary 
conditions. Unlike the expansion in t, we can in princi- 
ple calculate all the Z„ analytically, though this is not 
a terribly helpful procedure. What is important is that 



10 



the boundary conditions are trivial to implement - H44() 
will prescribe values for Z[^ at y = and y — \, which, 
together with Z{y,0) — Zq, fix an{t) and 6„(i) uniquely 
from (|7ni) . 

To show the existence of the T„, we need to examine 
the equation 



can be constructed as described previously, satisfying 



TvTt = 



L' 



Z'^n{Z) 



ZyZt 



2^V 



(71) 



which follows from the metric H45|l . T„ enters the LHS 
at 0(/i"), 

T,r,D/»"(T^ + S(y,t)) 

where again B{y,t) stands for binomial terms composed 
of already-known functions, the T,„ for to < n in this 
case. In fact the RHS of (|71|l does not contain Z„ at 
0{h"), since both Zy and Zt are 0{h). Therefore, the 
equation for r„ is of the form H59|) , specifically 

r;(y,t)=iC„(y,i) 

with the factor of t coming from the Zy. This can be 
integrated to 



T„(2/,t)=ytC„(2/,i)+d„(i) 



(72) 



where once again C„ is a binomial function of y and t. 
The arbitrary function dn{t) encodes the t reparameteri- 
sation invariance of the gauge. Because of the pre-factor 
y in ((7^ . the gauge T(0,i) = t corresponds to setting 
dn{t) = 0. We have therefore shown that there exist 
simple, binomial expressions for Tn{y,t) and Z„{y,t) for 
all n, with an extremely simple if tedious procedure for 
finding them explicitly. 

It should be emphasised, as in iJIVI that the precise 
nature of the boundary conditions are not important. 
Hence, although the explicit solution constructed in this 
section corresponds to empty, fine-tuned branes, it would 
be simple in principle to repeat the process for more gen- 
eral situations with both de-tuned branes and matter at 
the background level. 



VI. COSMOLOGICAL PERTURBATIONS 

We now consider more general spacetimes, where the 
spatial homogeneity is broken by cosmological perturba- 
tions. Our starting point is to consider a general back- 
ground spacetime, taking the metric to be of the form 

dsl = dY^ - n{Y, TfdT^ + a{Y, Tfn.jdx'dx^ (73) 

where flij is the relevant maximally-symmetric spatial 
metric of constant curvature. In this frame the branes 
are taken to have loci Y±{T). We then assume that a 
coordinate transformation 



9y9t 



i{f,g?fyft^O, 



(75) 



(76) 



with boundary conditions 

.g(y = 0,t)-y_(/(y = 0,t)), 

g{y^l,t)^Y+{f{y^l,t)), 

bringing the metric into the standard form 

dsl = d{tfdy^ - b{y, tfdt^ + a(y, tfVl.jdx'dx^ (77) 

with the branes fixed at y = and y = 1 (the precise 
forms of h(y, t) and a(y, t) are not important). Note that, 
in this section, we shall not fix the i-reparameterisation 
invariance as in (|21|l , which will allow us to consider more 
general situations than the branes moving apart after a 
collision, such as perturbations about static branes. 

We now allow the presence of cosmological perturba- 
tions, both to the background metric (|73|l and, in the case 
of scalar perturbations, to the positions of the branes 
themselves. We then attempt to find an infinitesimal 
gauge transformation which will not only retain the prop- 
erties 



55m = ^555 = 



(78) 



of the metric, but bring the brane positions back to y = 
and y = 1 if necessary. Firstly, we shall perform the coor- 
dinate transformation H74() to bring the background into 
the form H77I) , before constructing the gauge transforma- 
tion required to accommodate both the brane positions 
and the bulk metric perturbations into the SFB gauge at 
first order. Note that, in this section, we are still not con- 
sidering any non-gravitational degrees of the freedom in 
the bulk, which would in principle also need perturbing, 
such as a bulk scalar field. 



A. Tensor perturbations 

Firstly, we consider the simplest case of tensor pertur- 
bations. In five dimensions, the most general infinites- 
imal tensor perturbation [ig to H73|) only affects the 
transverse spatial components, i.e. 



9^j 



i{Y, T f {Vlij + hij) dx'dx' . 



(79) 



r = /(y,i), Y^g{y,t) 



(74) 



This will therefore not affect the properties (|78|) of the 
metric in the new frame (|77|l . Furthermore, the brane 
trajectories are scalars, and are therefore not affected by 
a purely tensor perturbation. Hence tensor perturbations 
can trivially be accommodated into the gauge H77|) . 
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B. Vector perturbations 

There are six vector degrees of freedom for five- 
dimensional perturbations about the background (|73|l . 
encoded in three divergenceless vectors: 



-,(T) 



(Y). 



ds^ = ds'o + Ql'dx'dT + QY'dx'dY (80) 
+a{Y,T)^F(^,yjdx'dx^ 

where | denotes the covariant derivative with respect to 
the metric flij and F,\ = Q)-. '^ = Q^ '^ = 0. Perform- 
ing the background coordinate transformation (|74|) . we 
obtain the perturbed version of (|77|l 



l(T). 



^(Y), 



ds' = ditydy'-biy,tydt' + [Q\'>fy + QY>gy)dyd 

+ (Qrft + QPgt)dtdx' 

+a{y,t)^n. 

Therefore, in order to satisfy the gauge conditions ()78|l . 
we require the hnear combination 



2i^(,|j-)) dxMx^ 



QPfv + QPgy 



(81) 



to vanish. Under the most general infinitesimal vector 
gauge transformation 

x'--^x'-+r]\ ?7|i = 0, 

the three vector perturbations transform as 

SQ'^P = -a{y,tfdyr,, 

5qP - -a{y,tfdtii^ 



SF, 



-m 



where rji — flij-rj^ . There are therefore two gauge- 
invariant vector perturbations (i.e. four degrees of free- 
dom) corresponding to Q^ — a^F- and Ql — a^Fi, and 
any one vector perturbation can be set to zero. In par- 
ticular, we can remove the linear combination (|81|l with 
suitable r]i satisfying 

fydyi], + gydtri^ = (qP fy + QP gy 

Again there is no vector perturbation to the brane tra- 
jectories, hence these still lie at y ~ Q and y = 1. 



Three of these can be removed at this stage by a scalar 
gauge transformation; effectively, this is just the state- 
ment that one can always choose a Gaussian-Normal 
gauge for the bulk where the off-diagonal terms Q\ , 
005 and the 1755 perturbation vanish: 



ds^ 



= ds, 



,-n'{Y,T)4>dT , ^_, 



Q^dTdx' 



^2a{y,tf{2i, 



E,,j)dx'dx^ 



(83) 



where (j) = 0oo- Performing the coordinate transforma- 
tion {75, we find 



ds^ = [d{t)'-n^f. 



{h^ + n^fh 



diN 



- 2n^fyft(l)dydt + Q^Jydx'dy 
Q^^ftdx'dt 



W [(1 + V')f^u + E^^\ dxMa;^' 



(84) 



where all the problematic terms have been written out 
in the first line, and n = n{y,t) = n{f{y,t),g{y,t)) etc. 
The background brane positions are fixed in this frame, 
but in the presence of scalar perturbations will be dis- 
placed to 

y+ (x, f) - 5y+ (x, t) , 2;_(x,f) = l + <5y_(x,t) (85) 

In five dimensions there are four gauge-invariant scalar 
degrees of freedom, hence we cannot eliminate any more 
of these perturbation variables with gauge transforma- 
tions. We can, however, perform another gauge trans- 
formation to rearrange them in such a way as to satisfy 
the conditions H78|) . whilst simultaneously bringing the 
two branes back to the fixed positions y — and y = 1. 
Under the most general infinitesimal scalar coordinate 
transformation 



the metric perturbations transform as 



(86) 



Sg55 = 


-2d^ 


i' -2ddx 


Sg5Q = 


-dH + b\' 


Sg5i = 


-d. 


a 



where ' — dy and ' — dt- Hence, from H84|l . we find the 
three conditions necessary to fix the required gauge: 



C. Scalar perturbations 

Scalar perturbations require a more detailed treat- 
ment, both because there are more of them in general 
and because one must also assume that the brane posi- 
tions are perturbed. In five dimensions there are seven 
scalar degrees of freedom. 



ds^ = dsl 



Assdr" + 2(/.n5dydr ~ n'iY, T)4>o^dT^ 



+Q^pdYdx' + Q^dTdx' 
+2a{y,tf{2i} + E,,j)dx'dx^ 



(82) 



2d^i" + 2ddx' +{11^ fl(l)) =0 

C + d^i - 2^C - Q/y = 
b\' ~ dH - n" fy ftcj) = 0, 



(87) 
(88) 
(89) 



together with the boundary conditions 

^(y = 0) = -<5y+(x,t), ^ (y = 1) = -5y_ (x, i) (90) 

which, from H85|) . will bring the brane positions to ?; = 
and y = 1. C and x have no boundary conditions to 
satisfy; once one has found the required ^, they follow 
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immediately from H88|l and (|89|l respectively. Eliminating 
x' from H87|) with 1)8 7|) , we obtain the equation of motion 
forC, 



e + k(y,t)i^F{y,t,^) 



(91) 



, , , dd 



2d2 



b^t 



fyft 



yJt' 



Interestingly only one of the original seven possible five- 
dimensional perturbations appears in this equation. In 
principle some information about the values of 4> and its 
derivatives at the boundaries can be inferred from the 
Israel junction conditions, though in practice this is not 
relevant. The perturbation is absorbed into a 'forcing' 
term F{y, t, x), with (|91(l taking the form of a generalised 
diffusion equation. 



1. Perturbations about static branes 

For static Randall-Sundrum branes (e.g. C — for 
empty branes), the original background configuration 
(|73|l is already in an SFB gauge with d{t) = 1. In this 
simple case, the equation of motion H91|l reduces to 



e"=F(y,i,x) 



(92) 



The general solution of this equation is 

^(y,i,x) = eiy,t,^) + (1 - 2/)A(x,t) +yB(x,t), (93) 

where S,^ is a particular solution to H92|l . and the linear 
term represents the most general homogeneous solution. 
A(x, t) and i?(x, t) are arbitrary functions which can be 
set so that ^ satisfies arbitrary boundary conditions at 
y = and y = 1. In particular, one can always solve the 
boundary conditions 1)90(1 with the choice 

^(x,t) = -(<52/+(x,i)+eP(0,i,x)) (94) 

S(x,t) == -(<52/-(x,t)+C^(l,t,x)) 

Hence any cosmological perturbations about a static, 
fixed-brane background can always be accommodated 
into an SFB gauge. 

2. Perturbations about general backgrounds 

For non-trivial k{y,t) the process of constructing a so- 
lution to (|91|l for arbitrary boundary conditions (|90f) is 
far more complicated and one must, once again, resort to 
a series solution. The boundary conditions (|90|l can be 
dealt with by a simple redefinition of the forcing function 
F{y,t,x), by setting 

e(2/,i,x) = w;(2/,i,x)-(l-y)<5+(x,t)-y^_(x,i). (95) 



The new function w has the simpler boundary conditions 
that w = at y = and y = 1, and satisfies the equation 

w" + k{y,t)w = G{y,t,x) (96) 

where 

G(y, t, x) = F{y, t, x)+fc(y, t) ((1 - y)6y+{x, t) + y 6y^{x, t)) 

We now expand the functions G, w and k in powers of 
t. Considering the case once again where the branes are 
moving apart after a collision which occurred at i = 0, 
we see from its definition in H91() that k{y, 0) = 0, hence 
there is no t" term in its expansion. We therefore write 

oo 

k{y,t) = 5]i"fc„(y) 

OC 

G{y,t,x) = ^i"G„(2/,x) 

oo 

w{y,t,x.) ^ ^rw„(y,x), 
t=o 

which, upon substitution into H96|) . gives 

w'^ - Go (97) 

ri-l 

w" 4- nkiWn = Gn- ^ mWm kn+l-m = G„ (98) 
m—1 

together with the boundary conditions 

w„(0,x) = u;„(l,x) =0. 

()98|l is just a simple, linear Sturm-Liouville problem 
which can always be solved with an eigenfunction ex- 
pansion subject to certain conditions. Specifically, we 
consider the eigenfunction equation 

C(2/) + A™fci(y)0™(y) = 0, 0„,(O) = 0„,(1) = 0. (99) 

The eigenfunctions 0m for m = 1, 2, 3... will form a com- 
plete, orthogonal basis, and so we can expand 

G„(y,i,x) = fci(j/)^G„™(x, i)0,„(y), 

m 
Wn{y,t,x) = y^W„m(x, t)0m(?/) 

771 

and these sums converge absolutely. Substituting these 
expressions back into 1)98(1 . we find 



Gr. 



W-n 



n- Xr. 



(100) 



and hence the required solution to 1(911) subject to the 
boundary conditions l(90() has been constructed. From 
((100(1 , we see that this solution only exists if the eigenval- 
ues of ((99(1 are either never integer- valued, or, if Am = N 



13 



for some integer iV, then the corresponding GNm van- 
ishes. 

Note that this conclusion is consistent with t- 
reparameterisation. For example, if we reparameterise 



t -^ T{t) = Ti f + T2 t^ 



then 



ki{y) 



kijy) 

N 



where N is the minimal positive integer such that T/y 7^ 
0. The eigenvalues of (|99|) will therefore all be multiplied 
by N, and Gnm is only non-zero for n a multiple of N. 
Hence the conclusion remains the same. This condition 
on the functions k{y,t) and G{y,t,:x.) can be regarded 
as necessary and sufficient for scalar perturbations to be 
accommodated within the SFB gauge. 

Note that if fci(y) = 0, the above argument will sim- 
plify considerably since each equation for Wn will just be 
of the form 

w'n — Gn for ki — 

which can be integrated immediately, with two arbitrary 
constants to fix the boundary conditions. Hence, in this 
case, it will always be possible to find the required gauge. 
However, for branes moving apart after a collision, from 
the above argument ki (y) will be proportional always to 
its value l/b{y, 0) in the gauge d{t) — t, which will not 
vanish. Hence one will have to check the eigenspectrum 
of the problem to verify the existence of the required SFB 
gauge. 

The precise nature of the function G{y, t, x) encoding 
both the bulk perturbation (j) and the perturbations Sy± 
to the brane positions is not important for the above dis- 
cussion. Hence one would expect, subject to the above 
caveats, to be able to include the effects of matter pertur- 
bations on the brane if desired. Since the perturbative 
construction of the background coordinate transforma- 
tion in ij iJIVI and Ivl were not sensitive to the precise na- 
ture of the boundary conditions either, one can expect 
to include matter at the background level as well. Hence 
we have presented in this paper a method for construct- 
ing SFB gauges very generally for Randall- Sundrum cos- 
mologies. 



VII. CONCLUSIONS 

The aim of this paper was to demonstrate the existence 
of gauges for Randall- Sundrum two-brane braneworlds in 
which the metric takes the particularly simple form ^ 
whilst, at the same time, the branes are located at fixed 
positions. Such gauges have been assumed to exist al- 
ready in the literature [22, |33, yj, Isll, without proof, 
since they simplify considerably the analysis of the five- 
dimensional Einstein equations necessary to derive ana- 
lytic results for effective theories. 



We first examined the dynamics of the system in the 
background, deriving from the trajectories of the branes 
the necessary boundary conditions for the coordinate 
transformation functions. These were then found using 
two different power series expansions, expanding either 
in time or in the parameter yC related to the mass of 
the bulk black hole. 

The former approach corresponds to an expansion in 
the radion, the proper distance between the two branes, 
and is adapted ideally to the case of two branes moving 
apart soon after a collision. This regime is of interest in 
recent attempts to model the B ig Bang as a collision be- 
tween two branes [11 113 111 111 Eliinil El . However, 
the non-linear equations which one needs to solve order 
by order in t do not admit an analytic solution. One 
can, however, demonstrate the existence of all terms in 
the series, satisfying the required boundary conditions, 
corresponding to being able to perform the necessary co- 
ordinate transformation correct to 0{t^) for arbitrary 
N. Whilst no analytic information about the conver- 
gence of the series is possible, numerical simulations show 
good agreement with the exact results, and we conjec- 
tured that the radius of convergence of the series would 
be equal to th, the time taken for the first brane to hit 
the horizon. 

A second method, namely an expansion in vC, cor- 
responds to an expansion in the velocity of the branes. 
It is possible to obtain analytic expressions for the coor- 
dinate transformation functions, and the resulting SFB 
metric, correct to arbitrary order in h. This is because 
the equations at each order are simple, linear differen- 
tial equations, involving just binomial functions which 
are trivial to integrate. 

Results on cosmological backgrounds are of limited in- 
terest; the important question is whether or not one can 
find an infinitesimal gauge transformation to accommo- 
date arbitrary cosmological perturbations, at first order, 
into an SFB gauge. We addressed this question last, and 
showed that this is trivially the case for tensor and vector 
perturbations. Scalar perturbations require considerably 
more work because one must also take into account a per- 
turbation of the brane trajectories themselves, and hence 
impose boundary conditions on the infinitesimal coordi- 
nate transformation ^ = 6x^. The second-order partial 
differential equation for ^ can be solved by an eigenvalue 
expansion, giving a well-defined constraint on the eigen- 
values of (|^ for the SFB gauge to exist. 

Whilst the arguments in this paper do not apply to 
braneworlds with a 'non-empty' bulk (i.e. with non- 
gravitational bulk degrees of freedom), they are not sen- 
sitive to the precise nature of the Israel junction condi- 
tions; one can include both detuned tensions and brane 
matter without altering the conclusions. 
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